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Abstract: 


1. A low storage scheme based on the universal matrices is developed for the discontinuous Galerkin time-domain method (DGTD). The 


proposed method can reduce the memory consumption dramatically with an acceptable raise in CPU time cost. 


2. The weighted Laguerre polynomials (WLP) scheme is integrated into the DGTD, thus leading to an unconditionally stable computational 
method. The resulted system uses the WLP as the temporal base and can be solved through many direct/iterative sparse linear system 
solvers. The WLP-DGTD method is suitable for the multi-scale simulations and eliminate the late-time instability of UPML in the 
DGTD method. 

3. 


Boundary Integral Method is combined with the WLP-DGTD method, and the computational accuracy versus mesh size and penalty 
factor 1s studied. 
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Numerical Flux 
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Semi-discrete form 
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Advantage Disadvantage 


e CFL constraint in explicit time 
discretization method 
* Large number of DOFs 














=) Challenges In DGTD 


Multiscale Problems ` Large Scale Problems 


Parallel Computation 


e Multi-core MPI-based Parallel Computation 

e Multi-core OpenMP-based Parallel 
Computation 

e Multi-GPU Acceleration Techniques (MPI- 
CUDA or MPI-OpenGL Hybrid Scheme) [10] 





Multiphysics and 
Non-linear 
Problems 
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Low Storage Scheme Based on Universal Matrices 
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"b Low Storage Scheme Based on Universal Matrices 
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Low Storage Scheme Based on Universal Matrices 
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Low Storage Scheme Based on Universal Matrices 


The comparison of memory cost and CPU 
time between low storage scheme and 
conventional scheme 
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`) Conclusion 


e We proposed a memory efficient scheme for DGTD based on 
universal matrix, and the memory usage is significantly reduced with 
slightly higher CPU time. 

* We integrated DGTD with Marching-On-in-Degree method. With this 


method, we eliminate the late-time instability of UPML in DGTD and 
overcome the CFL limitation imposed by the հ... 


* Boundary Integral Method is utilized in the DGTD-MOD method, 
numerical examples demonstrate the accuracy versus mesh size and 
penalty factor. 





ወ 


Tnank You! 





